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FLUX TRANSPORT SOLAR DYNAMO MODELS, CURRENT PROBLEMS 

AND POSSIBLE SOLUTIONS 
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RESUMEN 

Favor de proporcionar un resumen en espaiiol. If you are unable to translate your abstract into 
Spanish, the editors will do it for you. The sunspot solar cycle has been usually explained as the result 
of a dynamo process operating in the sun. This is a classical problem in Astrophysics that until the present is 
not fully solved. Here we discuss current problems and limitations with the solar dynamo modeling and their 
possible solutions using the kinematic dynamo model with the Babcock-Leighton approximation as a tool. In 
particular, we discuss the importance of the turbulent magnetic pumping versus the meridional flow circulation 
in the dynamo operation. 

ABSTRACT 

The sunspot solar cycle has been usually explained as the result of a dynamo process operating in the sun. 
This is a classical problem in Astrophysics that until the present is not fully solved. Here we discuss current 
problems and limitations with the solar dynamo modeling and their possible solutions using the kinematic 
dynamo model with the Babcock-Leighton approximation as a tool. In particular, we discuss the importance 
of the turbulent magnetic pumping versus the meridional flow circulation in the dynamo operation. 
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1. INTRODUCTION 

The sunspot cycle is one of the most interesting 
magnetic phenomenon in the Univ erse. It was dis- 



covered more than 150 years ago bv lSchwabd (|l844f ). 
but until now it remains an open problem in as- 
trophysics. There are several large scale observed 
phenomena that evidence that the solar cycle cor- 
responds to a dynamo process operating inside the 
sun. These can be summarized as follows: 

(1) The sunspots usually appear in pairs at both 
sides of the solar equator; the leading spot of a pair 
(i.e. the one that points to the E-W direction) has 
opposite polarity to the other one. Besides, leading 
spots in the northern hemisphere have the opposite 
polarity to that of the leading spots in the south- 
ern hemisphere. Sunspots invert their polarity ev- 
ery 11 years; the total period of the cycle is then 
22 years. This is known as the Hale's law; (2) A 
straight line connecting the leading and the com- 
panion spots of a pair has always an inclination of 
10° to 30° with respect to the equatorial line. This 
is known as the Joy's law; (3) When the toroidal 
field reaches its maximum, i.e., when the number of 
sunspots is maximum, the global poloidal field in- 



verts its polarity, so that there is a phase lag of 7r/2 
between the toroidal and poloidal components of the 
magnetic field; (4) The sunspots appear only in a 
belt of latitudes between ±30° (at both sides of the 
solar equator), these are known as the latitudes of 
activity; (5) The strength of the magnetic fields in 
the sunspots is around 10'^ G. The magnitude of the 
diffuse poloidal field is of tens of G. 



Parked (|l955fl was the first to try to explain the 
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solar cycle as a hydromagnetic phenomenon, since 
then although there has been important improve- 
ments in the observations, theory and simulations, 
a definitive model for the solar dynamo is still miss- 
ing. Helioseismology has mapped the solar internal 
rotation showing a detailed profile of the latitudi- 
nal and radial shear layers, which seems to confirm 
the usually accepted idea that the first part of the 
dynamo process is the transformation of an initial 
poloidal field into a toroidal one. This stage is known 
as the il effect. The second stage of the process, 
i.e.; the transformation of the toroidal field into a 
new poloidal field of opposite polarity is a less un- 
derstood process, and has been the subject of in- 
tense debate and research. Two main hypotheses 
have been formulated in order to explain the nature 
of this effect, usually denominated the a effect: the 
first one is based on the Parker's idea of a turbulent 
mechanism where the poloidal field results from cy- 
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clonic convective motions operating at small scales 
in the toroidal field. These small loops should re- 
connect to form a large scale dipolar field. How- 
ever, these models face an important problem: in 
the non-linear regime, i.e. when the back reaction 
of the toroidal field on the motions becomes impor- 
tant, the a effect can be cata strophically quenched 
(jVainshtein fc Cattaneo 1992) leading t o an ineffec- 
tive dynamo ( Cattaneo fc Hughe^ll996l PI. 

The second one is based on the formulations of 
Babcocd (|l96l[ ) and lLeightoi] (|l969[ ) (BL). They pro- 



posed that the inclination observed in the bipolar 
magnetic regions (BMR's) contains a net dipole mo- 
ment. The supergranular diffusion causes the drift 
of half of each of these active regions to the equa- 
tor and the drift of the other half in direction to the 
poles, so that this large scale poloidal structure an- 
nihilates the previous dipolar field. The new dipolar 
field is transported by the meridional circulations to 
the higher latitudes in order to form the observed po- 
lar field. This second mechanism has th e advantage 



of be i ng dir ectly observed at the surface (jWang et al 



19891 . [ 199l[ ). but it does not discard the existence of 



other a sources underneath. 

Following the BL idea, the physical model for the 
solar dynamo begins with a dipolar field. The differ- 
ential rotation stretches the poloidal lines and form 
a belt of toroidal field at some place within the so- 
lar interior, in the convective layer. This toroidal 
field is somehow pushed through the turbulent con- 
vective eddies and forms strong and well organized 
magnetic flux tubes. When the magnetic field is in- 
tense enough and the density inside a tube is lower 
than the density of the surrounding plasma, it be- 
comes unstable and begins to emerge towards the 
surface where it will form the BMRs. By diffusive de- 
cay, a BMR will form a net dipolar component with 
the opposite orientation of the original one, this new 
dipolar field is amplified during the cycle evolution 
until it reverses the previous dipolar field. In order to 
complete the cycle, it is necessary to transport this 
new poloidal flux first to the poles and then to the 
internal layers where the toroidal field will be cre- 
ated again, and so on. Most of the models in the BL 
mechanism use the meridional circulation flow as the 
main agent of transport. Recent works have invoked 



^Nevertheless, it is noteworthy that when the mag- 
netic heUcity is included in dynamical computations of o, 
it does not become catastrophycally quenched as long as 
the flux of the magnetic helicity remains non null (see 
[Brandenburg fc SubramanianI l2005al . for a complete review 
of this subject). The shear in the fluid could be the way 
through which the he licity flows to outside of the domain 
l lVishniac fc Choll200ll) . 



the turbulent pumping as an additional mechanism 
to advect the magnetic flux (see below). 

In the absence of direct observations to con- 
firm the model above, several numerical stud- 
ies have been performed in order to simu- 
late the solar dynamo. These can be divided 
in two main classes: global dynamical models 



els (iDikoati & Charbonneaul 
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first class integrates the full set of MHD equations 
in the solar convection zone and employ the inelastic 
approximation in order to overcome the numerical 
constrain imposed by fully compressible convection 
on the time-step. These models are able to reproduce 
the observed differential rotation pattern, but they 
do not generate a cyclic, and well organized pattern 
of toroidal magnetic field. The second class of sim- 
ulations solves the induction equation only and uses 
observed and/or estimated profiles for the velocity 
field and the diffusion terms. These models are rel- 
atively successful in reproducing the large scale fea- 
tures of the solar cycle, but the lack of the dynamical 
part of the problem has led to uncertainties in the 
dynamo mechanism. 

We have carried out several numerical tests with 
a mean field dynamo model in the Babcock-Leighton 
approximation in order to search for answers to four 
main issues: 

(1) where is the solar dynamo located?; (2) what 
is the dominant flux transport mechanism?; (3) how 
to explain the observed latitudes of the solar activ- 
ity?; and (4) why is the solar parity anti-symmetric? 

In the following paragraphs, we briefly summa- 
rize our model assumptions and, step by step, draw 
our approach to the questions above in the light of 
the numerical simulations. 

2. MODEL 

The equation that describes the temporal and 
spatial evolution of the magnetic field is the induc- 
tion equation: 

(9B 

'dt 

where U=Up + Jlr sin 6 is the observed velocity field, 
is the angular velocity, B=V x (Ae^) + B^e^, 
where Vx (Ae^) and are the poloidal and toroidal 
components of the magnetic field, respectively, 77 is 
the microscopic magnetic diffusivity and 



V X [U X B + 5 - r;V X B], 



1) 



e = aB + 7 X B - /3(V X B) 
- (5 X (V X B) - k(VB) , 



(2) 
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corresponds to the first and second order terms of 
the expansion of the electromotive force, u x b, and 
represents the action of the small-scale fluctuations 
over the large scales. The coefficients of ^ are 
the so-called dynamo coefficients. The first term 
on the right hand side of eq. (2) corresponds to 
the turbulent a effect coefficient, not considered 
in our Babcock-Leighton formulation. The second 
one is the turbulent magnetic pumping. The third 
corresponds to the turbulent diffusivity, which in 
our model is combined with the microscopic value 
{ijT—'H + P)- For the sake of simphcity, the other two 
terms are neglected. We solve equation ([T]) for A 
and Bif, with r and 9 coordinates in the spatial ranges 
O.6i?0— i?0 a nd 0— tt, respectively, in a 2 00 x 200 grid 
esolution fsee iGuerrero fc Munoa (2004), for details 
regarding the numerical model) . 

The profiles that we employ describe the re- 
sults of recent helioseismology inversions or numeri- 
cal simulations. For the differential rotation, we con- 
sider a profile mapped from helioseismology (see the 
continuous lines of Figure 1). For the meridional 
flow, we consider one cell per meridional quadrant, as 
usually assumed (doted lines in Figure 1). The alpha 
term (aB^) is concentrated between O.95i?0 and Rq 
and at the latitudes where the sunspots appear (see 
the continuous lines in Figure 2). Since it must result 
the emergence of magnetic flux tubes, we assume this 
term to be proportional to the toroidal field B^{rc, 0) 
at the overshoot interface between the radiative and 
the convective regions, rc=0.715i?o. For the mag- 
netic diffusion, we consider only one gradient of dif- 
fusivity located at rc which separates the radiative 
stable region (with ?7rz=10^ cm s~^) from the con- 
vective turbulent layer (with ?7cz=10^^ cm s~^) (see 
the dotted line in the upper panel of Figure 2). 

3. THE LOCATION OF THE SOLAR DYNAMO 
As remarked before, the differential rotation pat- 
tern which is responsible for the O effect, is revealed 
by high resolution helioseismology observations. It 
describes a solid-body rotation for the radiative core, 
and a differentially rotating convective layer with a 
retrograde velocity with respect to the radiative in- 
terior at higher latitudes and a pro-grade velocity 
at lower latitudes. The interface that bounds the 
solid-body rotation zone is named tachocline — its 
exact location and width have not been established 
yet. Another radial shear layer has been recently 
identified just below th e solar photosphere in the up- 
per 35 Mm of the sun (jCorbard fc Thompson|[200a ) 
(see the gray dashed line in Fig. 1). With this 
newly discovered shear layer, it is even more diffi- 
cult to define where the dynamo operates. There 



o 




Fig. 1. Isorotation lines of tlie solar angular velocity 
as inferred from helioseismology observations (continu- 
ous lines) and the adopted meridional circulation stream 
lines. The red dot-dashed line shows the boundary be- 
tween the overshoot layer and the convection zone. 



has been so far, an apparent common agreement that 
the dynamo is operating at the tachoc line. Howev er 



this possibility has several problems (jBrandenburs 
20051) . One of the main difficulties is that toroidal 
flux ropes formed in the tachocline should have in- 
tensities ~ lO^'-lO'' G in order to become buoy- 
antly unstable and to emerge at the surface to 
form a BMR the appropria t e tilt g iven by th e Joy's 
law (iD'Silva fc Choudhuril Il993t iFan et all [l993t 



Caligari et al.ril995l ll998riFan fc Fisheil Il99a iFan 
20041 ). One important limitation of this scenario is 
that 10^ G results an energy density that is an or- 
der of magnitude larger than the equipartition value, 
so that a stable layer is required to store and am- 
plify this magnetic field. This raises another ques- 
tion with regard to the way in which the magnetic 
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Fig. 2. Radial and latitudinal profiles for a (continu- 
ous line), 77t (dotted line) and for the pumping terms 7r 
and 7£) (dashed and dot-dashed lines, respectively). All 
the profiles are normalized to their maximum value. Ex- 
tracted from \ Guerrero & de Gouveia Dal Pind \200^) . 

flux is dragged down to deeper layers. In r ecent work 
(jGuerrero fc de Gouveia Dal Pind l2007a[ ). we have 
explored the contributions of the shear terms in the 
dynamo equation, i^p-V)Q.=Brd^/dr + B0/rd^d9 
with the aim of determining where the most strong 
toroidal magnetic fields are produced. We found 
that the radial shear component is about two or- 
ders of magnitude smaller than the latitudinal com- 
ponent. Therefore, when a new toroidal field begins 
to develop its growth is dominated by the latitudi- 
nal shear. Its amplification begins in the bulk of the 
convection zone and it is transported to the stable 
layer where it must reach the desired magnitude (see 
Fig. 3). How this fiux is transported is the subject 
of the next section. 

If a near-surface shear layer is turned on in the 
model (Guerrero & de Gouveia Dal Pino 2008 and 
references therein) , two main branches appear in the 
butterfly diagram (Fig. 4). One is migrating pole- 
ward (at high latitudes) and another is migrating 
equatorward (below 45°). This resu lt is expected 
if the Parker- Y oshimura sign rule ( Parker 1955t 
Yoshimura Il975l ) applies. This contribution to the 




Fig. 3. (a) Butterfly diagram and latitudinal snapshots 
for the toroidal (b) and the poloidal (c) fields. The dark 
(blue) and light (red) gray (color) scales represent pos- 
itive and negative toroidal fields, respectively; the con- 
tinuous and dashed lines represent the positive and neg- 
ative poloidal fields. Only toroidal fields greater than 
2 X 10* G (the most external contours) are shown in 
panels (a) and (b). This model started with an anti- 
sym metric initial condition (see §6 for d e tails) . Extracted 
from \Guerrero & de Gouveia Dal Pind \200^} . 



toroidal field was prev iously explored in a B abcock- 
Leighton dynamo by iDikpati et al. (2002). They 
have discarded the near-surface radial shear layer 
because it generates butterfly diagrams in which a 
positive toroidal field gives rise to a negative radial 
field, which is exactly the opposite to the observed. 
Our results, on the other hand, present the correct 
phase lag between the fields. This difference prob- 
ably arises from the fact that wc arc using a lower 
meridional circulation amplitude. Anyway, as can be 
seen in Fig. 4, the polar branches are strong enough 
to generate also undesirable sunspots close to the 
poles. The period increases to 15.6 y due to the fact 
that the dominant dynamo action at the surface goes 
in the opposite direction to the meridional flow (i.e, 
the dynamo wave direction is dominating over the 
meridional flow). 
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Fig. 4. Butterfly diagram for a model with the same pa- 
rameters as in Fig. 3, but with near-surface shear layer. 
For this model T=15.6 yr B^^^ Jr=0.715) = l.l x 10^ 
5.^.™aJr=0.98)=1.9 X 10" G and Br„„, =131.7 G. This 
model start ed with symmetric initial conditi o ns. Ex- 
tracted from \ Guerrero & de Gouveia Dal Find \200A ). 

4. FLUX TRANSPORT MECHANISMS 

F or an ctO dynamo, the Parker - Yoshimura sign 
rule (Parker 1955t F^shimura 1975 ) establishes that 
the direction of the dynamo wave is equatorward or 
poleward if the product a ■ dil/dr is < or >0, respec- 
tively. Hence, models with a solar like rotation law 
operating at the tachocline, with a positive a effect 
in the northern hemisphere, as believed, and without 
meridional circulation, should result a solution with 
magnetic bran ches migrating in t he opposite way to 



that observed (jKiiker et al.ll200lh . Models with one 



cell of meridional circulation, poleward at the sur- 
face and equatorward at the base of the convection 
zone, produce result s with the appropriate direc- 
tion of propagatio n ( Dikpati fc Charbonneau 19991 : 



Kiiker et all l200ll : iBonanno et all 120021 ). It was 



demonstrated also that t he meridional circulation 
sets t he period of the cycle ( Dikpati fc Charbonneaul 
19991) and that it is the most logical way to trans- 
port the novel poloidal fields to the inner layers 
in the dynamo process. These models in which 
the time of the cycle fits better with the advec- 
tive time than with the diffusive one are usually 
called advection dominated or flux-transport dy- 
namos. There is, however, an important problem 
with these models: they require a large scale merid- 
ional flow and this is observed only at the surface. 
It is possible that the real meridional flow is too 
weak at the inner regions to p enetrate the over- 
shoot layer and the tachocline (jCilman fc Miesch 
2004; Riidiger et al.l 2005), or perhaps, it has a 



The turbulent pumping effect corresponds to the 
transport, in all directions, of magnetic flux due to 
the presence of density (buoyancy) and turbulence 
(diamagnetism) gradients in convectively unstable 
layers. In the FOSA (First Order Smoothing Ap- 
proxim ation, see e.g. [Brandenburg fc Subramanian 
l2005b . and references therein), the radial component 
of the diamagnetic pumping can be calculated as- 
suming a linear dependence with the variations in 
the magnetic diffusivity Udia=—'^VT/'2, the buoy- 
ancy component depends on the density gradients 
Ubuo=—ff^p/p- In the boundary between the solar 
convection zone and the overshoot layer, it is prob- 
able that the diamagnetic velocity is of the order of 
50 ms-i (jKitchatinov fc Riidigeill2008f) . This value 



strengthens the importance of the pumping relative 
to the assumed radial meridional flow velocity which 
is < 10 ms^^. For our model, we obtain J7dia=— 47 
cm s~^ when a variation of two orders of magni- 
tude is cons idered in the diffusivity in a thin region 
of 0.015i?o (jCuerrero fc de Gouveia Dal Pincllioosl ) 
(the buoyancy component of the pumping is not con- 
sidered because the density is not a parameter in 
kinematic models). 

The effects of turbulent pumping have been 
rarely considered in mean field dynamo models. A 
first approach showing its importance in t he so - 



lar cycle was made by [Brandenburg et al.l (jl992i ): 



since then few works have incorporated the dia- 
magnetic pumping component in the dynamo equa- 
tion as an extra diffusive term that c an provide a 
downward vel ocity, as discussed above ( Kiiker et al.l 



2001 



Bonanno et all [2002. 20Qi). More recently. 



Kapvla et al.l ( 2006b[ ) have implemented simula- 
tions of the mean field dynamo in the distributed 
regime, including all the dynamo coefficients pre- 
viousl y evaluated in rnagneto-convection simula - 
tions ( Ossendriiver et al. 20021 : Kapvla et al.ll2006a ). 



They produced butterfly diagrams that resemble the 
observations. However, to our knowledge no special 
efforts have been made to study the pumping effects 
in the meridional plane (i.e., inside the convection 
zone) or in a BL description. 

We have included the turbulent pumping terms 
(see the dashed and dot-dashed lines in Fig. 2 
which correspond to the radial, 7^, and latitu- 
dinal, 7e, pumping terms) calculated from local 



multicell pattern (iMitra-Kraev fc ThompsonI l2007t magneto-convection simul ations (|Ossendrijver et al 



Bonanno et"al]|2006l : Ijouve fc Brunll2007t l. If this is 
the case, it is very hard to explain the equatorward 
migration of the toroidal branches and it is neces- 
sary to find another flux-transport mechanism. The 
turbulent pumping seems to be a good candidate. 



20021 iKapvla et al.ll2006al) into the induction equa- 
tion (eq. 1). 

For a dynamo model operating at the tachocline, 
we find that the pumping terms lead to a distinct lat- 
itudinal distribution of the toroidal fields when com- 
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Fig. 5. Meridional flow streamlines and the butterfly 
diagram for a model with the full pumping term, but 
with a shallow meridional flow penetration with a depth 
of only O.8i?0, [/o=1300 cm s~^, 7eo=90 cm s~^ and 
7ro=30 cm s~^. For this model we obtain r=10.8 yr, 
S0„.„=4.5 X fO* G and Br„„, =154.9 G. This model 
start ed with anti-symmetric initial conditions . Extracted 
from \ Guerrero & de Gouveia Dal Pind If200i }. 

pared with the results of Fig. 3. The turbulent and 
density gradient levels present in the convectively un- 
stable layer cause the pumping of the magnetic field 
both down and equartorward, allowing its amplifica- 
tion within the stable layer and its later emergence 
at latitudes very near the equator. 

If we also include in this model recent helioseis- 
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Fig. 6. The same as in Fig. 5, but for a model 
with near-surface shear action. For this model T—16.3 
yr ■B^„oJr=0.715)=9.7 x 10'', Jr=0.98)=1.9 x 

10" G and -B^^^^ =164.4 G. This model started 
with symmetric initial conditions. Extracted from 
\Guerrero & de Gouveia Dal Pindl \200A) . 



mic results ([Mitra-Kraev fc Thompsonl 120071 ) that 
suggest that the return point of the meridional circu- 
lation can be at ^ O.95-R0, at lower regions, beneath 
~ O.8-R0, a second weaker convection cell or even a 
null large scale meridional flow can exist. In Figure 
5, we obtain a butterfly diagram that agrees with the 
main features of the solar cycle, besides, wc find that 
in this case, it is the pumping terms that regulate 
the period of the cycle, leading to a different class 
of dynamo that is advection-dominated by turbulent 
pumping rather than by a deep meridional flow. 

For the dynamo model operating at the near- 
surface shear layer, we find that the toroidal fields 
created at high latitudes are efficiently pushed down 
before reaching a significant amplitude, so that only 
the equatorial branches below 45° survive (see Fig. 
6). This result is explained by the fact that the ra- 
dial pumping component has its maximum ampli- 
tude close to the poles (see the dashed line in Fig. 
2). These results also agree with the observations. 
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5. HOW TO EXPLAIN THE OBSERVED 
LATITUDES OF THE SOLAR ACTIVITY? 

The radial shear at the tachochne has its maxi- 
mum amphtude in regions close to the poles, for this 
reason, it is a common problem in mean field dynamo 
models to present large unde sirable toroidal mag- 



9x1 0* 



netic fi elds in the polar regions. iNandv fc Choudhuri 



(120021 ) proposed a deep meridional flow as a way to 
avoid the formation of strong toroidal fields at high 
latitudes. Under this assumption, the strong toroidal 
magnetic fields formed at high latitudes are pushed 
down, inside the radiative zone, by the meridional 
flow and are stored there until they reach latitudes 
below 30° . However this assumption may lead to un- 
desirable mixing of the chemical elements between 
the radiative and convective zones and to problems 
regarding the angular momentum transfer. Besides, 
some results of numerical simulations show that the 
meridional flow is unable to penetrate the tac hocline 
(|Gilman fc Mieschll2004HRiidiger et al.ll2005h . 

We have explored this subject in two ways. 
First, we built a hybrid mod el in which we com 
bined the profiles used bv iNandv fc Choudhuri 



(2002) with those used bv iDikpati fc Charbonneaij 

19991) and allowe d a deep meridional flow 

Guerrero fc Mufiozl 12004 ). We found that the 
high-latitude toroidal fleld is sensitive to the 
model. Then, we explored this problem by 
changing the sh ape and the thickness of the 



solar tachocline ( Guerrero fc de Gouveia Dal Pinol 
20073^ 1 and have found that the thinner the 
tachocline, the smaller the intensity of the toroidal 
magnetic field at high latitudes (see Fig. 7). A 
thin tachocline must be fully contained inside the 
overshoot zone, in such a way that only part of 
the poloidal magnetic field is able to reach it and 
then produce a small quantity of toroidal field. The 
toroidal field generated there is not strong enough 
to emerge. In Fig. 8, we compare the toroidal fields 
produced for a thin tachocline with those produced 
with an intermediate one. This result is also depen- 
dent of the magnetic diffusivity in the convection 
zone, as it can be seen in Fig. 7, however, for a thin 
tachocline the models always result weak toroidal 
fields above 60°. 

With regard to this question of the distribu- 
tion of the toroidal field, the role of the pumping 
is also important. While it pushes poloidal fields 
inside the tachocline, it also pumps all the field 
equatorward, in such a way that the toroidal fields 
formed inside the convection zone due to latitudi- 
nal shear will go fast to the stable region where 
they are stored and amplified before the eruption 
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Fig. 7. Maximum of the toroidal magnetic field at the 
top of the tacholine as a function of the diffusivity (in 
log-scale) at a latitude of 60°. The different line styles 
represent different widths of the tachocline di . The dot- 
ted line represents the limit between buoyant and non- 
buoyant magnetic fields 5 x 10* G. Only the values below 
this line will appear at the desired la t itudes Extracted 
from \Guerrero & de Gouveia Dal Find t2007d) . 



(jCuerrero fc de Gouveia Dal Pinoll200^ . 

In the case that we consider the sunspots as the 
product of toroidal fields being formed at the near- 
surface layer, the latitude of activity is easily ex- 
plained by the pumping, as described in the previous 
section (see Fig. 6). 

6. THE PARITY PROBLEM 

The anti-symmetric parity observed in the so- 
lar cycle is one of the most challenging questions 
in the solar dynamo theory. The magnetic parity 
in a model may depend on the location of the a ef - 
fect ( Dikpati fc Gilmaiill200ll ; iBonanno et al. I l2002h . 
or on the diffusive couplin g between the poloida l 



2004^ 



field in both hemispheres (jChatteriee et al 
but it may also be the result of the imprint 



of 



the quadrupolar form of the meridion al flow on the 
poloid al magnetic fleld, as argued bv ICharbonneaul 
(j2007l ). This could explain why models with the a ef- 
fect located at the upper layers (where the magnetic 
Reynols number and the quadrupolar imprint of the 
meridional flow are larger) tend to a quadrupolar 
parity faster than models with the a effect located 
at the tachocline. 

We have explored a little this problem in order to 
see whether the turbulent pumping plays some role 
in it. We have found that models without pump- 




Fig. 8. Snapshots of the positive (negative) toroidal field contours are shown in blue (red) scale together with the field 
hnes of the positive (negative) poloidal field in continuous (dashed) lines, for four different times (T/8, r/4, ST/S and 
r/2) along a half 11-year cycle. The upper (a) and middle (b) panels present the same snapshots for models with a 
thin (O.O2_R0) and an intermediate (O.O6-R0) tachocline Extracted from \Guerrero & de Go uveia Dal Pirw (2007i). 



ing result a quadrupolar solution. When beginning 
with a dipolar initial condition, they spend several 
thousand years before switching to a quadrupolar 
solution (see Fig. 9a ). This result div e rges f rom 
the one obtaine d by 



Dikpati fc GilmanI (|2001l ) or 



IChatteriee et al. ( 2004 ) in which the change begins 
only after around 500 yr. This result indicates the 
strong sensitivity of the parity to the initial condi- 
tions, in such a way that, for example, the present 
parity observed in the sun could be temporary. 

The models with full pumping (e.g.. Fig. 5) 
conserve the initial parity either it is symmetric or 
anti-symmetric (see Fig 9b). This suggests that the 
strong quadrupolar imprint due to meridional cir- 
culation can be washed out when the full turbulent 
pumping is switched on. 

On the other hand, the full pumping in models 
that include near-surface shear tend to the dipolar 
parity from the first years of integration (see con- 
tinuous line of Fig. 9c). We explain this result as 
a product of a better coupling between the poloidal 
fields in both hemispheres (this coupling is due to 
the local a effect considered in these cases), plus the 
action of the pumping eliminating the effect of the 
quadrupolar component of the meridional flow on 
the poloidal magnetic field. 



7. SUMMARY AND DISCUSSION 

We have used a mean field dynamo model in the 
BL approach in order to look for answers to four 
current problems widely reported in the literature of 
solar dynamo modeling. 

Our results confirm the idea that there should be 
a magnetic layer below the convection zone where 
magnetic fields which are mainly produced in the 
convection zone are stored. We have found that it is 
possible to have a flux-transport dynamo without a 
well defined meridional flow pattern, dominated by 
the pumping advection. However, efforts in order 
to obtain a more realistic profile for the meridional 
velocity profile arc still required, as well as to obtain 
a better comprehension of the real contribution of 
the pumping velocities. 

Nevertheless, numerical simulations including 
near-surface shear as observed, provide also support 
for a near-surface magnetic layer, since toroidal fields 
with intensities between 10^ - 10^ G are formed there 
and since the radial differential rotation is negative 
at lower latitudes, the direction of migration of the 
butterfly wings constructed with these flelds repro- 
duces the one observed (e.g.. Fig. 6). 

The latitude of activity established by the ob- 
servations (between ±30°) can be explained by both 
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a) Without pumping 
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time steps number (33 kyr) 

t)) With full turbulent pumping 
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time steps number (44 kyr) 

c) With near — surtace shear 
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Fig. 9. Parity curves for the three classes of models con- 
sidered, i.e., (a) for models without pumping (as, e.g., 
in Fig. 3); (b) for models with full pumping (as, e.g., 
in Fig. 5); and (c) for models with near-surface shear 
(as, e.g., in Figs. 4 and 6). In the panels (a) and (b), 
the continuous, dashed and dot-dashed lines correspond 
to symmetric, anti-symmetric and random initial condi- 
tions, respectively. In panel (c), the continuous line is 
used for the model with turbulent pumping while the 
dash ed line is for the model without pumping Extracted 
f row ] Guerrero & de Gouveia Dal Pind \200t ). 



scenarios above for the magnetic layer. For magnetic 
fields being stored in the overshoot region, a thin 
(< O.2i?0) tachocline could be the solution to avoid 
strong toroidal fields at the polar regions. On the 
other hand, if the near-surface shear is considered, 
then the magnetic pumping provides the required 
downwards flux of the weak polar fields, letting only 
the toroidal fields to survive at the active latitudes. 

The parity in a dynamo solution is a problem that 
requires especial attention. We have investigated 
this problem looking for the role that the pump- 
ing could play in the solutions. Our simulations 
support the idea that the quadrupolar solution in 
most of the models is due to the strong quadrupo- 
lar imprint due to the one-cell meridional fiow pat- 
tern. This imprint is larger at the surface than at 
the bottom of the convection zone, therefore the 



result s of iDikpati fc GilmanI (|200ll ) ; iBonanno et al 
(j2002i ) which suggest that an a effect operating at 
the tachocline results an anti-symmetric solution, 
could be explained by this fact. We suggest that 
models with an a effect operating close to the sur- 
face are also able to generate anti-symmetric solu- 
tions, as observed, if a mechanism such as pumping 
cleans the quadrupolar imprint. 

More observational and theoretical efforts are 
necessary in order to determine where the feet of 
the magnetic flux tubes responsible for the sunspots 
are located. This is an issue that needs to be solved 
before a more realistic coherent dynamo model can 
be constructed. 
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